Magnetic and transport percolation in diluted magnetic semiconductors 
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The ferromagnetic transition in a diluted magnetic semiconductor with localized charge carriers is 
inevitably a percolation transition. In this work we theoretically study the correlation between this 
magnetic percolation and transport properties of the sample, including the possibility of a simultane- 
ous transport percolation. We find nontrivial signatures of the percolating magnetic clusters in the 
transport properties of the system, including interesting non-monotonic temperature dependence of 
the system resistivity. 
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I. INTRODUCTION 

"Diluted magnetic semiconductors" (DMS) is a com- 
mon name for a wide class of semiconductor materials 
doped with magnetic impurities (with their relative con- 
centration up to a few percent) which become ferromag- 
netic at low temperatures. Magnetic semiconductors are 
now actively investigated both theoretically and experi- 
mentally, due to their potential applications in new gen- 
erations of semiconductor devices. 

An example of a DMS material is given by 
Gai_a;Mn2;As (GaAs doped with Mn). For the ferromag- 
netic transition to occur, the concentration of Mn impuri- 
ties should be relatively high, x ~ 0.01 — 0.05. Such high 
concentration of impurities makes growth of magnetic 
Gai_a;Mn2:As a difficult problem, which was solved only 
relatively recentlyi— In order to prevent phase separation, 
magnetic Gai_a;Mna;As should be grown at low temper- 
atures (T ~ 200 — 300°C), which results in abundance 
of different types of crystal defects. These problems are 
not specific to Gai_a;Mna;As, in fact, all presently known 
magnetic semiconductors suffer from the same problems. 
As a result, theoretical study of magnetic semiconduc- 
tors is very difficult: two factors - strong disorder and 
exchange interaction (which causes ferromagnetic transi- 
tion) - must be taken into account non-perturbatively. 

The mechanism of ferromagnetic transition is common 
for most (if not for all) magnetic semiconductor systems. 
Since the concentration of magnetic impurities is only of 
the order of few percent, their direct exchange interaction 
does not play the decisive role in determining the ground 
state of the system. The ferromagnetic transition occurs 
due to interaction between impurity spins mediated by 
charge carriers, which are present in appreciable quanti- 
ties because of the high concentration of defects in the 
crystal lattice. This picture seems to be accepted more or 
less universally, but the detailed nature of the transition 
varies greatly from a system to a system. Because of the 
large number of different magnetic semiconductor ma- 
terials, high defect concentration in them, and possible 
variations of growth processes, it is highly unlikely that a 
"universal theory of magnetic semiconductors" will ever 
be developed. A theory limited to a certain subclass of 



magnetic semiconductor materials, which share same ba- 
sic properties, is more likely to provide useful information 
on the nature of the processes underlying the ferromag- 
netic transition and give some insight into experimental 
results. 

In this paper, we study ferromagnetic semiconductors 
that in the limit of zero temperature become insulating. 
The mechanism of ferromagnetic transition for an ex- 
treme case of such a system, when the charge carriers are 
pinned to some point defects, was considered in our ear- 
lier paperii^ Here we include transport properties in our 
consideration, and also extend our scope to the systems 
where the charge carries are localized within extended re- 
gions rather than pinned with the wave function decaying 
exponentially from the localization center. 

A characteristic feature of DMS physics is the pres- 
ence of a metal-insulator transition as a function of the 
magnetic impurity concentration x. Typically for low x 
{< 0.03) the system (e.g. Gai_j,Mnj,As) is an insula- 
tor with strongly localized carriers and with the resis- 
tivity increasing exponentially as T ^ 0. For larger x 
the system tends to be "weakly" metallic in the sense 
that the resistivity, although very large (~ 0.01 fl ■ cm), 
does not exhibit strong temperature dependence. There 
are reports^ of a reentrant insulating phase at higher Mn 
concentration, x > 0.05 (More recent experiments^ on 
carefully annealed samples do not find this reentrant in- 
sulating behavior at higher Mn concentration.) In all 
situations both metallic and insulating phases are found 
to be ferromagnetic. Other DMS systems show typically 
the same behavior (i.e. the existence of a metal-insulator 
transition with a weakly metallic state, which has a very 
high resistivity, for higher values of a;) or often just an 
insulating ferromagnetic state without any metallic be- 
havior whatsoever. Systems showing a metal-insulator 
transition tend to have higher ferromagnetic transition 
temperature Tc in their metallic phases, and therefore 
understanding the insulating ferromagnetic DMS phase 
takes on significance, not just from the perspective of 
fundamental physics involving interplay of disorder and 
magnetism, but also from the practical technological mo- 
tivation of enhancing Tc in DMS materials. This is par- 
ticularly true since the current zeroth-order theoretical 
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understanding of the optimally doped metallic DMS fer- 
romagnetism is based on a simple continuum virtual crys- 
tal approximation Weiss mean-field theory treatment of 
an RKKY-Zener model^ that is intrinsically invalid in 
the insulating regime. In fact, it is currently unclear 
whether the observed DMS metal-insulator transition is 
a disorder-driven impurity band or valence band tran- 
sition. It is, however, clear that disorder plays a criti- 
cal role in DMS physics, both in its transport properties 
and in the magnetic properties. We focus in this work 
on understanding DMS transport properties in the lo- 
calized insulating regime by applying our recently devel- 
oped theory)^ which has been successful in explaining the 
magnetic properties on the insulating side of the metal- 
insulator transition. We hope to elucidate the role of dis- 
order and magnetic clustering on DMS transport prop- 
erties through this theoretical study. Although the DMS 
systems with the highest Tc do not belong to this class, 
their resistivity is still very large, so disorder must play 
a very important role in the properties of these systems. 
Study of DMS systems that are insulating at T — > 0, 
as done in this paper, can shed some light on this role, 
which is difficult to understand using the approach based 
on the picture of free charge carriers. 

So the systems we consider in this paper are insulating 
at T = 0, which means that the Fermi level is below the 
mobility edge in the impurity band or below the bottom 
of the conduction band (for electrons; for holes it would 
be the valence band) if there are no extended states in the 
impurity band. From now on, we will refer to this class 
of magnetic semiconductor systems as "systems with lo- 
calized charge carriers." Taking into account the high 
defect concentration in a typical magnetic semiconduc- 
tor material, we must conclude that the charge carrier 
density in systems with localized charge carriers is highly 
inhomogeneous. Since it is charge carriers that transmit 
exchange interaction between magnetic impurities, this 
effective interaction must also be highly inhomogeneous. 
Thus, when the temperature is lowered, the ferromag- 
netic transition will first occur locally, in the regions with 
higher charge-carrier density, i.e. with stronger effective 
exchange interaction between magnetic impurities. As 
the temperature goes further down, these finite-size re- 
gions, which have random sizes and positions, will grow 
and merge, until finally the ferromagnetic correlation is 
established across the whole sample.^ Such a scenario, 
which must hold for any magnetic semiconductor with 
localized carriers, implies that ferromagnetic transition is 
inevitably a percolation transition, with clusters of fer- 
romagnetic regions already existing at T > Tc and an 
infinite magnetic cluster opening up at T = Tc- This 
percolation picture of ferromagnetism in insulating DMS 
systems is now well-accepted and has been verified-'^ via 
direct numerical simulations. 

The foundation for the study of ferromagnetic percola- 
tion transition was laid down in the 70s for diluted ferro- 
magnetic alloys, where magnetic atoms interact directly 
with each other. ^ Early numerical simulations of local- 



ized DMS systems^^ already indicated indirect evidence 
for magnetic percolation playing a role in the magnetic 
transition. Our earlier papeB^ extended results of Ref. 
to magnetic semiconductor systems, where the exchange 
interaction between magnetic atoms is not direct, but is 
rather mediated by charge carriers. 

The word "percolation" by itself assumes some trans- 
port. Even though the mathematical percolation the- 
ory with its physical applications goes beyond this lit- 
eral meaning (an example is given by diluted ferro- 
magnetic alloys of Ref. 9, where no transport is con- 
sidered), one may still naturally wonder if a percola- 
tion transition in a given system means some enhance- 
ment of charge transport. In ferromagnetic semiconduc- 
tors, this question is quite relevant, since charge carri- 
ers, which have some ability to move around the sam- 
ple, are present, and since the ferromagnetic percolation 
transition is facilitated by these very carriers. In fact, 
some experiments on magnetic semiconductors did ob- 
serve connections between ferromagnetic transition and 
enhancement of charge transport, up to the point of 
the temperature dependence of resistivity being non- 
monotonic with the maximum around the ferromagnetic 
transition temperature In addition, the existence 
of the metal-insulator transition in DMS systems makes 
transport percolation considerations relevant since there 
is an intrinsic connection between metal-insulator transi- 
tion and percolation. We emphasize in this context that 
the magnetic percolation transition, considered in Ref.|3, 
is purely a statement on the percolation properties of the 
temperature-dependent magnetic clusters in DMS mate- 
rials, and has little to do (in a direct sense) with transport 
percolation properties (which must take into account the 
site-to-site hopping of the localized charge carriers, as we 
do in this paper) - it is, in principle, entirely possible 
for a system to have a magnetic percolation transition 
with no transport percolation altogether. In fact, these 
are two separate physical phenomena whose relationship 
is theoretically explored in this paper in the specific con- 
text of DMS systems. 

The main purpose of this work is to establish con- 
nection between ferromagnetic transition in a magnetic 
semiconductor system with localized carriers (which, as 
we mentioned above, must have percolation nature) and 
transport properties of such a system. The outline of 
the paper is as follows: In Sec. ^ we introduce our 
model. In Sec. IIIII we consider the case of charge car- 
riers strongly pinned to localization centers. The ferro- 
magnetic properties in such a system were studied in our 
earlier papers(Sii^ here we discuss its transport proper- 
ties. The considerations of this section are based on the 
impurity positions being uncorrelated. The case of impu- 
rities forming clusters in considered in Sec. II VI Finally, in 
Sec. we qualitatively discuss the transport properties of 
a system, whose charge carriers are on the verge of being 
delocalized. For all these systems, we also discuss depen- 
dence of the resistivity on the applied magnetic field. 
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II. MODEL HAMILTONIAN 

Charge carriers in Gai-^^Mnj^As are holes, donated by 
Mn impurities. While other magnetic semiconductors 
may have electrons as the carriers mediating exchange 
interaction between magnetic impurities, in this paper 
we will refer to the charge carries as "holes," just for the 
sake of brevity. All our conclusions would hold no matter 
what the charge carriers are. 

In fact, the whole model system we use is based on 
Gai-rcMuajAs, where Mn atoms act both as magnetic 
impurities and acceptors, thus providing both local mo- 
ments, which order ferromagnetically at T < Tc, and 
charge carriers, which mediate ferromagnetic interaction 
between these local moments. We do not expect our 
results to be specific to this particular model as far as 
the qualitative results are concerned, so our conclusions 
should be valid for other systems carrier-mediated insu- 
lating DMS as weU. 

In our model, a Mn impurity in Gai_a;Mna;As is pre- 
sented by two spin-degenerate levels, which we will refer 
to as "the deep level" and "the shallow level." The deep 
level, when occupied by one electron, provides the impu- 
rity's spin. Coulomb repulsion prevents a second electron 
from entering this level. The shallow level plays the role 
of the impurity's acceptor level. Since the electron wave 
function of this level has larger localization radius, the 
Coulomb interaction between a deep-level electron and a 
shallow-level one is weaker than that between two elec- 
trons on the deep level. 

The full Hamiltonian of the system reads: 



^0 = E 



a 



mn^mct^na. 



mnoc 

V (t^^'^'^c^ a +hc 



(1) 



The last term of Eq. describes hopping from the 
deep level of one impurity to the shallow level of an- 
other. Hopping between deep levels is neglected because 
of the rapid fall-off of the electron wave functions at these 
levels. We take the Fermi energy equal to zero, £i^- = 
throughout the paper. 

The parameters of Hamiltonian ^ must be chosen to 
mimic magnetic impurities in Gai_3;Mn2;As. The deep 
level in the ground state must be taken by only one hole, 
so we take Ud > £d > 0. The shallow, acceptor levels 
of the impurities are to donate, on average, less than 
one hole per impurity, so the Fermi level coincides with 
the energy of a hole placed onto the shallow level of an 
impurity, whose lower level level is already taken by one 
hole, which means that Uod—^o = £f = 0. The bottom of 
the conduction band is assumed to be separated from the 
impurity level Uod — ^o by the energy Eq 3> T, tmn, so the 
(dclocalized) valence band states can be safely excluded 
from our consideration. 

We assume the wave functions of these localized lev- 
els to fall off exponentially, with the characteristic de- 
cay radii equal to ad and ao for deep and shallow levels 
respectively, with ad oq. Then the hopping matrix 
elements tmri and i„„, which are proportional to the 
overlap of the corresponding wave functions, are given 
by 



ao / \ ao ' 



ao 



(2b) 



With these parameters, we may make the Schrieffer- 
Wolff transformationi^ for the deep levels, thus reducing 
Hamiltonian to 



in 

+ E* 



Uc)(^ -t-C -^C^ I C I 
^ m 7 m t m [ m i 



mn ma na 



lmna.(3 



(3) 



where a = [ux, o'yjCrz) is the vector of Pauli matrices. 



where indices m and n run over magnetic impurities, aj^^^ 
and cj^^ are the creation operators for an electron with 
spin a localized at the mth impurity at the deep/shallow 
level respectively. The first and the third terms in the 
brackets in Eq. represent the energies of holes occu- 
pying the deep and the shallow level respectively. The 
second and the fourth terms in the brackets describe 
the Coulomb repulsion between two holes on each of 
these two levels, while the fifth term is for the inter-level 
Coulomb repulsion. The first term after the brackets ac- 
counts for hole hopping between the shallow levels of two 
impurities. 



ji _ AOd) (od) / J_ 1 

•^mn ^ ''7711 ''In \ ^ ^ tt 



and 



exp I )(4) 



ai/3 



(5) 



The other parameters of the system are the concentration 
of impurities rij and the concentration of holes rih. The 
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spatial distribution of impurities is assumed to be random 
throughout the paper, although not necessarily uncorre- 
lated. We also assume, as it is the case in Gai_a;Mnj;As 
and many other DMS, that the system is very heavily 
compensated, rih rii, with many more magnetic impu- 
rities that charge carriers present in the system. 

Depending on the values of the system parameters, the 
system described by Hamiltonian ijjjl may be either in- 
sulating or metallic in the limit T ^ 0. In this paper we 
will concentrate on the former case. Ferromagnetic tran- 
sition in such a system inevitably has percolation nature, 
as we pointed out in the introduction. 

III. STRONGLY LOCALIZED CHARGE 
CARRIERS 

In this section, we consider the case of aQUi 1, with 
the impurity positions being uncorrelated. Randomness 
in the impurity positions leads to randomness in the hop- 
ping matrix elements tmn, which exponentially depend on 
the distances between the corresponding impurities. 
This strong (by orders of magnitude, provided agni <IC 1) 
variation of tmn creates strong variation of localization 
potentials for the holes, which thus will be tied to the 
pairs of least-separated impurities. We must note that 
the concrete mechanism of hole localization does not have 
crucial impact on the derivations and conclusions pre- 
sented in this section. In order to keep our presentation 
straightforward we limit ourselves to the framework of 
Hamiltonian ©, when the holes are localized at these 
closest impurity pairs. Any additional disorder can be 
taken into account by inclusion of new terms to Eq. (jS)) 
and proper modification of the hole localization parame- 
ters. 

With strongly localized holes, for study of magnetic 
properties of the system Hamiltonian ||2J) can be reduced 
to: 

H = ^ ^ JmjSmSji (6) 

mj 

where indices m and j label magnetic impurities and 
holes respectively, and Sm/^j are the impurity/hole spin 
operators. The expression for the matrix elements J^j 
of the impurity-hole exchange interaction depend on the 
hole localization mechanism. In the case of localization 
by the pairs of least-separated impurities, it is given by 

Jmi = 7: { J C-) ^ J (2) (7) 

2 V / 
for a hole localized at the impurities labelled by indices 
n^"'^^ and ^T-j^^ 

Since the concentration of holes is much smaller than 
that of impurities, rih ^ ni, one localized hole is sur- 
rounded by many magnetic impurities. Exchange inter- 
action between the (localized) holes and magnetic impu- 
rities leads to their mutual polarization when tempera- 
ture T is below exchange constant Jmj- Since exchange 



coupling Jmj decays with the impurity- hole distance rmj , 
see Eq. (j2b|l and the first impurities to get their 
spins aligned with a hole's spin as the temperature de- 
creases are the ones most close to the hole's localization 
site. At lower temperatures, more distant impurities have 
their spins polarized by the hole to whose domain they 
belong. This complex consisting of a hole and magnetic 
impurities polarized by it is conventionally called "bound 
magnetic polaron." The characteristic radius size of a 
polaron grows as the temperature decreases!^ 

rpo\ (T) = Y In ^7-^ , (8) 

where A ^ 1 and 

^-('rr(^-i^). 

cf. Eq. As polarons overlap, they form polaron 

clusters with all impurities belonging to a given cluster 
having their spins aligned in the same direction. The 
temperature^i 

Tc - sS\Jo\ (ajjnh) V'^iM exp ( ^h^] (g) 

V (ai]"h)V 

at which the infinite cluster spanning the whole sample 
appears is the ferromagnetic transition temperature. The 
quantitative theory of this transition was recently devel- 
oped and presented by us in Ref. Qj and further developed 
m Ref. hj. In this paper we concentrate on transport 
properties of such a system and their correlation with 
the magnetic properties. 

At temperatures of the order of the ferromagnetic tran- 
sition temperature, the hole transport in the system oc- 
curs by means of nearest-neighbor hopping. The resis- 
tivity p of such a system depends on temperature expo- 
nentially, 

p a exp(£;hop/r) . (10) 

The characteristic hopping activation energy ii^hop has 
two contributions. First, there is random level mismatch 
between two localization sites, produced by the disorder. 
The second contribution to the hopping activation en- 
ergy comes from the interaction of a localized hole with 
neighboring magnetic impurities. This interaction leads 
to polarization of impurity spins by the hole spin and low- 
ers the energy of the resulting bound magnetic polaron. 
The energy associated with this polarization (polaron's 
"binding energy") is given by 

£^poi = STrJoao'^i , (11) 

provided rpoi > oq, which is satisfied at temperatures of 
the order and below the ferromagnetic transition temper- 
ature Tc- When a hole hops to another localization site, 
it abandons the region it has polarized and lands among 
non-polarized impurities unless this new localization site 
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is within another bound magnetic polaron. The "binding 
energy" Epd of a polaron is much less than the charac- 
teristic energy level mismatch E^^J^ ~ to at different lo- 
calization sites. Therefore this "binding energy" will not 
have a noticeable effect on the hopping trajectory of a 
hole. At temperatures close to the ferromagnetic transi- 
tion temperature Tc, while the bound magnetic polarons 
barely touch each other, it means that most of the hops 
will be to unpolarized regions, see Fig. Therefore, the 
characteristic hopping activation energy will acquire an 
addition term Spoi, 



hop 



= E, 



(0) 
hop 



^poi at T < 



(12) 



At low enough temperatures, however, the infinite clus- 
ter of bound magnetic polarons will cover the whole sam- 
ple, so wherever a hole hops, it will land in a region which 
is already polarized in a optimal way, so the impurity 
polarization will not have any contribution to the hop- 
ping activation energy. Our Monte-Carlo studies show 
that the characteristic temperature Tcovcr at which the 
infinite cluster covers most of the sample corresponds to 

^pol(-^covcr) ^ 

2rpoi(T'c), which yields 



TcCxp -- 



0.86 



Thus 



E, 



hop 



<l at T < r„ 



(13) 



The concrete scenario of the infinite cluster spread- 
ing over the transport paths highly depends on the sam- 
ple, because even though these two random networks are 
not completely independent, they generally do not coin- 
cide. Since the size of a bound magnetic polaron grows 
logarithmically slowly as the temperature goes down, 
sample-dependent variations in the mutual arrangement 
of the transport paths and the infinite cluster lead to 
strong variations in the temperatures at which the tran- 
sition from the high-temperature [Eq. (|12|l ] to the low- 
temperature [Eq. H13(l ] exponent occurs and in the pre- 
cise functional form of p(T) during this transition. How- 
ever, since the transition is slow, the overall temperature 
dependence p{T) of the resistivity must be monotonic, 
because even though magnetization reduces the hopping 
barrier, the decrease in temperature make it harder for 
electrons to get activated to overcome it. 

External magnetic field B applied to the sample polar- 
izes the impurities all over the sample, with the average 
value of the impurity polarization given by 



(5,) = SBs 



where /is is the Bohr magneton 
Bs{x) 



2s + 1 , 2s + 1 
coth — X 



2s 



2s 



1 , 1 
— coth — X 
2s 2s 



(14) 



(15) 




FIG. 1: Charge transport in a DMS system with strongly 
localized carriers. The spins of localized holes are shown with 
black arrows, the regions with magnetic impurities polarized 
by these spins are dotted. Beyond the percolation thresh- 
old, the polarized regions coalesce, forming an infinite cluster 
spanning the whole sample. The charge transport is facil- 
itated by holes hopping between the localization sites, (a) 
On-site repulsion does not allow a hole to hop onto a site al- 
ready taken by another hole, (fe) A hole can hop only to a free 
site, which may be outside the region polarized by other holes, 
(c) Hops to already polarized regions are rare at T ~ Tc. 




FIG. 2: Sketch of the temperature dependence of resistivity 
in a DMS system with strongly localized carriers. The hop- 
ping activation energy i?hop is given by Eqs. 1121 and 11311 for 
regions (1) and (2) respectively, with the transition between 
the regions being sample-dependent. Even though Shop is 
lower is region (2), the temperature decrease needed to re- 
duce it is large enough for p(T) to remain monotonic. 



is the Brillouin function, and the direction of the z axis is 
chosen to be parallel to the magnetic field. As the result 
of this polarization, the impurities around the localiza- 
tion site a hole hops into are not entirely uncorrelated 
with those around the site it hops from. The "polaron 
binding energy" Hll|) is therefore reduced by 



AEp,i{B) = -STTJoalmBs 



gif^BB \ 



(16) 



and the resulting magnetic-field dependence of the resis- 
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tivity is given by 
p{B) = p\b=oCxp 



Fermi level is 



ksT I 



(17) 



as it follows from Eqs. ifTUIl . ((T^ . and ((TE|) . One can 

see that the resistivity decreases as the applied magnetic 
filed grows, in agreement with the experimental results^ 



IV. CHARGE CARRIERS LOCALIZED WITHIN 
CLUSTERS 

The previous section addressed transport in magnetic 
semiconductors with strongly localized charge carriers. 
"Strongly localized" in that case meant that a charge 
carrier is pinned to some point (a pair of close impuri- 
ties or some other defect in crystal structure), with the 
carrier localization radius being smaller than the char- 
acteristic distance between pinning centers. Coulomb re- 
pulsion prevents a carrier from entering a localization site 
already taken by another carrier. This regime takes place 
when concentrations rij of magnetic impurities and rih of 
charge carriers are small. For such small n; and rih, fer- 
romagnetic transition temperature is typically of the 
order of few Kelvins for realistic parameters. For practi- 
cal applications, however, it is desirable to have Tc higher 
than few Kelvins, which implies higher ni and Tih, so one 
may naturally wonder what happens to the charge trans- 
port when the magnetic impurity concentration is so high 
that the assumptions of Sec. lIIII are no longer applicable. 

If the positions of impurities are uncorrelated, charge 
carriers become delocalized when parameter rijag exceeds 
some critical value of the order of unity. The analyti- 
cal description of this transition is very challenging, and 
hardly any quantitative results can be obtained. We 
postpone the discussion of this case until the next sec- 
tion. In this section we consider the case when magnetic 
impurities are grouped into clusters. This arrangement 
is not improbable at all, taking into account significant 
differences between the non-magnetic atoms of the host 
lattice and the magnetic dopants. In fact there are the- 
oretical indicationgi^ that the impurities should group 
into clusters, due to electrostatic interactions during the 
growth process. 

For the parameter region in the vicinity of nia^g ^ 1 
clustering of impurities means that a hole will be able to 
move freely within a cluster, but not between clusters, i.e. 
the characteristic size of a hole wave function will be of 
the order of a cluster size rather than of the size of a wave 
function of a hole localized at an isolated impurity or an 
impurity pair. Because of the hole wave function being 
spread over many impurities, the Coulomb energy which 
repels a hole from a site occupied by another hole is re- 
duced as compared to the case of strongly localized holes 
considered in Sec. IIIII If the cluster size is large enough, 
several holes may enter one cluster, and the characteris- 
tic magnitude of the Coulomb energy in a cluster at the 



£^coui ^ U 



(18) 



where A^h and iVj are the number of impurities in the 
cluster and the number of holes sitting in it, respectively. 
Since rii ^ rih, for a typical cluster we must have Ni 3> 
Nh, and, therefore i?coui ^ U- 

The activation energy i?hop for hole hopping between 
the clusters contains contributions from the Coulomb- 
and disorder-induced level mismatch between the clus- 
ters and from the exchange interaction between the holes 
hopping between the clusters and the impurities forming 
these clusters. The first contribution stems mainly from 
the Coulomb repulsion between the holes in clusters. The 
characteristic value of the gap between the highest occu- 
pied and the lowest unoccupied levels in a cluster is given 
by Eq. H18|) . The relative positions of these levels in dif- 
ferent clusters are completely random, so a characteristic 
value of the mismatch between the highest occupied and 
lowest unoccupied levels in two neighboring clusters will 
also be of the order of i?coui given by Eq. H18|) . The ran- 
dom mismatch Egp between two single-particle levels in 
two clusters also plays some role, but the characteristic 
value of this mismatch 



E. 



tn 



sp 



^2/3 1/3 



(19) 



should be much less than -Ecoui for realistic systems. 

The second contribution to the hopping activation en- 
ergy comes, as it was stated above, from the exchange 
interaction between holes and impurities. At high tem- 
peratures, spins of impurities and holes in a cluster are 
not polarized, and there is no contribution to the hopping 
activation energy due to exchange. At lower tempera- 
tures, however, exchange interaction between magnetic 
impurities and holes leads to ferromagnetic transition. 
The mechanics of the latter is as follows: first, exchange 
interaction between holes and impurities in clusters they 
occupy lead to ferromagnetic transition within clusters. 
The exchange interaction between the clusters is expo- 
nentially small, oc exp(— Tci/ao), where rd is the charac- 
teristic distance between the clusters. Therefore, as the 
temperature is lowered and ferromagnetic transition oc- 
curs within clusters, the magnetic moments of different 
clusters are still uncorrelated. At even lower tempera- 
tures, the correlation between the clusters is established, 
and the macroscopic ferromagnetic transition occurs. At 
high temperatures when the even clusters are not polar- 
ized a hole hopping from a cluster to another one expe- 
riences no change of the exchange energy. However at 
temperatures close to the ferromagnetic transition tem- 
perature the clusters are inevitably polarized, and this 
creates an additional energy barrier for a hole hopping 
between two uncorrelated clusters. 

While the ferromagnetic transition in a cluster does 
not imply full polarization of all impurities in it, at the 
temperatures of the order of the ferromagnetic transition 
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FIG. 3: Charge transport in a DMS system with clustered 
magnetic impurities. The picture refers to the temperatures 
low enough for the ferromagnetic transition to occur within 
clusters, whose spins are shown with black arrows. Holes hop 
between the clusters preferring shorter gaps, so the percola- 
tion paths coincide with transport paths. 




T 

FIG. 4: Sketch of the temperature dependence of resistivity 
in a DMS system with clustered impurities. Since percolation 
and transport paths coincide, ferromagnetic transition leads 
to immediate decrease in the hopping energy, which may re- 
sult in non-monotonic temperature behavior of the sample 
resistivity. 

temperature all the impurities must be polarized, and the 
energy associated with the exchange interaction of holes 
in a cluster with all the impurities in it equals 

Ecxch — Jo ■ (20) 

When two clusters are not correlated, a hole hopping 
from one cluster to the other one experiences the change 
of potential of the order of i^cxch, so the expression for 
the characteristic hopping activation energy reads 

Ehop — -Ecoui + -Ecxch at Tc < T . (21) 

At temperatures below Tc, the correlation between the 
cluster spins is established, and a hole hopping into an- 
other cluster experiences exactly the same exchange po- 
tential as in the cluster it has left. This reduction of the 



exchange contribution to the hopping activation energy 
from Jo to zero (roughly speaking) occurs while tem- 
perature T is still of the order of Tc- Therefore, the 
exponent in Eq. H1U|I may actually become smaller as 
the temperature goes down past provided the second 
term in Eq. (|21|l . which vanishes when the spins of clus- 
ters become aligned, is larger than the first, temperature- 
independent term, that is 

< Jo (22) 

The external magnetic field applied to a sample low- 
ers its resistance, similarly to the case discussed in the 
previous section. However, now the reduction of the re- 
sistivity in the systems is due to polarization of impurity 
clusters as a whole, not of individual impurities, since 
all impurities in a cluster have their spins aligned in the 
same direction. Since polarization of a cluster by mag- 
netic field of a given magnitude strongly depends on the 
size of the cluster, and since the sizes of impurity clusters 
are random with possibly wide distribution, the quanti- 
tative dependence of resistivity on the applied magnetic 
field is highly sample-dependent. 

V. CHARGE CARRIERS ON THE VERGE OF 
DELOCALIZATION 

At sufficiently high concentrations of impurities and 
holes, some of the holes are delocalized even in the limit 
T — > 0. The ferromagnetic transition in systems of this 
class does not have percolation nature, unless the frac- 
tion of delocalized carriers is very small. With physics of 
their ferromagnetic transition and charge transport being 
entirely different from that of the systems with localized 
charge carriers, these systems are beyond the scope of 
the current paper. In this section, we consider systems 
which are close to localization threshold, but still not be- 
yond it. In the ground state of such a system, all charge 
carriers are localized within some finite regions, and the 
resistivity of the sample goes to infinity as T ^ 0. How- 
ever, because of the small separation between the Fermi 
level and the mobility edge, finite temperature excites 
some holes into delocalized states. It is hardly possible 
to make any quantitative statements about magnetiza- 
tion and resistivity of the system, but still some qualita- 
tive statements can be made with the help of intuition 
we have developed dealing with the case of the previous 
section. 

Since, despite some holes being thermally excited into 
delocalized states, many holes are still localized, the ef- 
fective exchange interaction between the magnetic impu- 
rities, which is induced by holes, is still highly inhomoge- 
neous across the sample. There will be "puddles" filled 
with holes, where interimpurity exchange concentration 
is strong, and relatively hole-free regions, which are cov- 
ered only by tails of hole wave functions decaying away 
from these puddles, see Fig.|Sl Qualitatively, this picture 
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is similar to that of Sec. lIVI shown in Fig.|31 with the dif- 
ference that in the case considered in the present section 
transport paths go mostly along hole-filled regions with 
fewer gaps on the way. These gaps are bridged by the 
delocalized states, so there is no need for hopping. 

Similarly to the case of previous sections, it is re- 
gions with higher hole density that undergo (local) fer- 
romagnetic transition first, with ferromagnetic correla- 
tions across the sample established at lower tempera- 
tures. When two neighboring "hole puddles" have the 
different spin orientation, a delocalized hole going from 
one puddle to the other encounters an additional barrier, 
which can partially reflect it back or even entirely pre- 
vent it from entering the puddle. When the spins of two 
puddles oriented in the same direction, as it is the case 
at the temperatures below the ferromagnetic transition 
temperature, no such problem arise, and a delocalized 
hole move freely from a puddle to a puddle across the 
sample. 

We thus have two competing phenomena that occur 
as the temperature goes down. On one hand, at lower 
temperatures we have fewer holes excited to delocalized 
states. On the other hand ferromagnetic transition oc- 
curing as temperature goes down makes it easier for the 
hole to travel across the sample. The interplay of these 
two phenomena determines the temperature dependence 
of resistivity. Due to the extremely complex nature of 
hole wave functions close to the mobility edge, we are un- 
able to provide any qualitative results. However, with the 
system under consideration being in some sense an ex- 
treme case of the system consider in the previous section, 
we may deduce that the non-monotonic behavior of resis- 
tivity must be possible, and that the non-monotonicity 
becomes more pronounced as the Fermi level of the holes 
approaches the mobility edge, i.e. when the concentra- 
tions of holes and impurities become higher. Similarly to 
the previous cases, magnetic field applied to the system 
would align the spins of hole puddles and enhance the 
transport. 

VI. CONCLUSION 

We have considered correlations between ferromag- 
netic percolation transition and charge transport in mag- 
netic semiconductors with localized charge carriers. We 
found that ferromagnetic transition alone would enhance 
charge transport, but the temperature decrease, which 
is needed for the transition, may completely eliminate 
these gains. In a system with strongly localized charge 
carriers (Sec. decrease in the hopping rate of the 

charge carriers due to temperature decrease overcomes 
decrease in the hopping activation energy due to ferro- 
magnetic transition, and the temperature dependence of 
the resistivity is always monotonic. At higher concen- 
trations of magnetic impurities and charge carriers, the 
situation may be the opposite, and resistivity may ex- 



perience a dip as the temperature goes down past the 




FIG. 5: Charge transport in a DMS system with charge 
carriers on the verge of being delocalized. Dotted areas show 
"puddles" of charge carriers, with impurity spins correlated 
across a puddle even at temperatures above the ferromagnetic 
transition. Gaps between puddles are bridged by thermally 
excited delocalized states (hatched region) . 



ferromagnetic transition temperature (Sees. HVIIVl) . The 
higher the concentrations of the impurities and carriers 
are, the more pronounced the dip is. In all cases consid- 
ered in the paper, the resistivity of the system decreases 
if an external magnetic field is applied. 

It is important to point out that our theoretical DMS 
transport results are in qualitative agreement with the 
available experimental transport results in the localized 
insulating regime, where the resistivity increases strongly 
(exponentially) with temperature as T — > 0. In par- 
ticular, deep into the insulating regime, far away from 
the metal-insulator transition point, the experimentally 
measuredSiii resistivity rises monotonically with decreas- 
ing temperature similar to our results (Fig. ^ presented 
in Sec. IIIII for strongly localized charge carrier system. 
On the other hand, close to the metal-insulator transition 
(but still on the insulating side) transport experiments in- 
deed observa^iii a striking temperature-dependent non- 
monotonicity in the measured resistivity as we find in 
our theory for "not-so-strongly" localized charge carriers 
(SecHVl Fig. EI). Unfortunately, the insulating DMS fer- 
romagnetic regime has not yet been extensively studied 
experimentally (compared with the metallic DMS regime 
with optimal Tc values) although our view is that this 
regime is as important and as interesting as the metallic 
regime in terms of the development of our understand- 
ing of DMS physics. We urge that more experimental 
work done in the insulating ferromagnetic regime (and 
on the metal-insulator transition itself) for us to develop 
a better and more quantitative understanding of trans- 
port and magnetization phenomena in DMS systems. 
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